use "C:\Users\andre\Desktop\Artigos\Artigo_Ens_Sup_2\Bases\Agregadas_Stata\Agregada.dta", clear

estimates drop _all


********************************Modelo para conclusão do EM

*Sem controles
logistic ens_medio b(first).ano##b(7).cat_chefe_agg if (rec_mod==1) & ((ano==1992 | ano==2002 | ano==2022) | (ano==2012 & versao==1)) [pweight = peso] 

*Com controles

*1992
quietly logistic ens_medio b(7).cat_chefe_agg b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==1992) & (rec_mod==1) & (((ano>=1992 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_cat_1992

quietly margins cat_chefe_agg, post
estimates store margins_ens_medio_cat_1992

quietly nlcom (CAT1_EM: _b[1.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT2_EM: _b[2.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT3_EM: _b[3.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT4_EM: _b[4.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT5_EM: _b[5.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT6_EM: _b[6.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT8_EM: _b[8.cat_chefe_agg] / _b[7.cat_chefe_agg]), post
estimates store nlcom_ens_medio_cat_1992


*1993
quietly logistic ens_medio b(7).cat_chefe_agg b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==1993) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_cat_1993

quietly margins cat_chefe_agg, post
estimates store margins_ens_medio_cat_1993

quietly nlcom (CAT1_EM: _b[1.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT2_EM: _b[2.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT3_EM: _b[3.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT4_EM: _b[4.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT5_EM: _b[5.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT6_EM: _b[6.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT8_EM: _b[8.cat_chefe_agg] / _b[7.cat_chefe_agg]), post
estimates store nlcom_ens_medio_cat_1993

*1995
quietly logistic ens_medio b(7).cat_chefe_agg b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==1995) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_cat_1995

quietly margins cat_chefe_agg, post
estimates store margins_ens_medio_cat_1995

quietly nlcom (CAT1_EM: _b[1.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT2_EM: _b[2.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT3_EM: _b[3.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT4_EM: _b[4.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT5_EM: _b[5.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT6_EM: _b[6.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT8_EM: _b[8.cat_chefe_agg] / _b[7.cat_chefe_agg]), post
estimates store nlcom_ens_medio_cat_1995

*1996
quietly logistic ens_medio b(7).cat_chefe_agg b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==1996) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_cat_1996

quietly margins cat_chefe_agg, post
estimates store margins_ens_medio_cat_1996

quietly nlcom (CAT1_EM: _b[1.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT2_EM: _b[2.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT3_EM: _b[3.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT4_EM: _b[4.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT5_EM: _b[5.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT6_EM: _b[6.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT8_EM: _b[8.cat_chefe_agg] / _b[7.cat_chefe_agg]), post
estimates store nlcom_ens_medio_cat_1996

*1997
quietly logistic ens_medio b(7).cat_chefe_agg b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==1997) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_cat_1997

quietly margins cat_chefe_agg, post
estimates store margins_ens_medio_cat_1997

quietly nlcom (CAT1_EM: _b[1.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT2_EM: _b[2.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT3_EM: _b[3.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT4_EM: _b[4.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT5_EM: _b[5.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT6_EM: _b[6.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT8_EM: _b[8.cat_chefe_agg] / _b[7.cat_chefe_agg]), post
estimates store nlcom_ens_medio_cat_1997

*1998
quietly logistic ens_medio b(7).cat_chefe_agg b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==1998) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_cat_1998

quietly margins cat_chefe_agg, post
estimates store margins_ens_medio_cat_1998

quietly nlcom (CAT1_EM: _b[1.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT2_EM: _b[2.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT3_EM: _b[3.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT4_EM: _b[4.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT5_EM: _b[5.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT6_EM: _b[6.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT8_EM: _b[8.cat_chefe_agg] / _b[7.cat_chefe_agg]), post
estimates store nlcom_ens_medio_cat_1998

*1999
quietly logistic ens_medio b(7).cat_chefe_agg b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==1999) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_cat_1999

quietly margins cat_chefe_agg, post
estimates store margins_ens_medio_cat_1999

quietly nlcom (CAT1_EM: _b[1.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT2_EM: _b[2.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT3_EM: _b[3.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT4_EM: _b[4.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT5_EM: _b[5.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT6_EM: _b[6.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT8_EM: _b[8.cat_chefe_agg] / _b[7.cat_chefe_agg]), post
estimates store nlcom_ens_medio_cat_1999

*2001
quietly logistic ens_medio b(7).cat_chefe_agg b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2001) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_cat_2001

quietly margins cat_chefe_agg, post
estimates store margins_ens_medio_cat_2001

quietly nlcom (CAT1_EM: _b[1.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT2_EM: _b[2.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT3_EM: _b[3.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT4_EM: _b[4.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT5_EM: _b[5.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT6_EM: _b[6.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT8_EM: _b[8.cat_chefe_agg] / _b[7.cat_chefe_agg]), post
estimates store nlcom_ens_medio_cat_2001

*2002
quietly logistic ens_medio b(7).cat_chefe_agg b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2002) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_cat_2002

quietly margins cat_chefe_agg, post
estimates store margins_ens_medio_cat_2002

quietly nlcom (CAT1_EM: _b[1.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT2_EM: _b[2.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT3_EM: _b[3.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT4_EM: _b[4.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT5_EM: _b[5.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT6_EM: _b[6.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT8_EM: _b[8.cat_chefe_agg] / _b[7.cat_chefe_agg]), post
estimates store nlcom_ens_medio_cat_2002

*2003
quietly logistic ens_medio b(7).cat_chefe_agg b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2003) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_cat_2003

quietly margins cat_chefe_agg, post
estimates store margins_ens_medio_cat_2003

quietly nlcom (CAT1_EM: _b[1.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT2_EM: _b[2.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT3_EM: _b[3.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT4_EM: _b[4.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT5_EM: _b[5.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT6_EM: _b[6.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT8_EM: _b[8.cat_chefe_agg] / _b[7.cat_chefe_agg]), post
estimates store nlcom_ens_medio_cat_2003

*2004
quietly logistic ens_medio b(7).cat_chefe_agg b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2004) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_cat_2004

quietly margins cat_chefe_agg, post
estimates store margins_ens_medio_cat_2004

quietly nlcom (CAT1_EM: _b[1.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT2_EM: _b[2.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT3_EM: _b[3.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT4_EM: _b[4.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT5_EM: _b[5.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT6_EM: _b[6.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT8_EM: _b[8.cat_chefe_agg] / _b[7.cat_chefe_agg]), post
estimates store nlcom_ens_medio_cat_2004

*2005
quietly logistic ens_medio b(7).cat_chefe_agg b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2005) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_cat_2005

quietly margins cat_chefe_agg, post
estimates store margins_ens_medio_cat_2005

quietly nlcom (CAT1_EM: _b[1.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT2_EM: _b[2.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT3_EM: _b[3.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT4_EM: _b[4.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT5_EM: _b[5.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT6_EM: _b[6.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT8_EM: _b[8.cat_chefe_agg] / _b[7.cat_chefe_agg]), post
estimates store nlcom_ens_medio_cat_2005

*2006
quietly logistic ens_medio b(7).cat_chefe_agg b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2006) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_cat_2006

quietly margins cat_chefe_agg, post
estimates store margins_ens_medio_cat_2006

quietly nlcom (CAT1_EM: _b[1.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT2_EM: _b[2.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT3_EM: _b[3.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT4_EM: _b[4.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT5_EM: _b[5.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT6_EM: _b[6.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT8_EM: _b[8.cat_chefe_agg] / _b[7.cat_chefe_agg]), post
estimates store nlcom_ens_medio_cat_2006

*2007
quietly logistic ens_medio b(7).cat_chefe_agg b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2007) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_cat_2007

quietly margins cat_chefe_agg, post
estimates store margins_ens_medio_cat_2007

quietly nlcom (CAT1_EM: _b[1.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT2_EM: _b[2.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT3_EM: _b[3.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT4_EM: _b[4.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT5_EM: _b[5.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT6_EM: _b[6.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT8_EM: _b[8.cat_chefe_agg] / _b[7.cat_chefe_agg]), post
estimates store nlcom_ens_medio_cat_2007

*2008
quietly logistic ens_medio b(7).cat_chefe_agg b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2008) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_cat_2008

quietly margins cat_chefe_agg, post
estimates store margins_ens_medio_cat_2008

quietly nlcom (CAT1_EM: _b[1.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT2_EM: _b[2.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT3_EM: _b[3.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT4_EM: _b[4.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT5_EM: _b[5.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT6_EM: _b[6.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT8_EM: _b[8.cat_chefe_agg] / _b[7.cat_chefe_agg]), post
estimates store nlcom_ens_medio_cat_2008

*2009
quietly logistic ens_medio b(7).cat_chefe_agg b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2009) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_cat_2009

quietly margins cat_chefe_agg, post
estimates store margins_ens_medio_cat_2009

quietly nlcom (CAT1_EM: _b[1.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT2_EM: _b[2.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT3_EM: _b[3.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT4_EM: _b[4.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT5_EM: _b[5.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT6_EM: _b[6.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT8_EM: _b[8.cat_chefe_agg] / _b[7.cat_chefe_agg]), post
estimates store nlcom_ens_medio_cat_2009

*2011
quietly logistic ens_medio b(7).cat_chefe_agg b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2011) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_cat_2011

quietly margins cat_chefe_agg, post
estimates store margins_ens_medio_cat_2011

quietly nlcom (CAT1_EM: _b[1.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT2_EM: _b[2.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT3_EM: _b[3.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT4_EM: _b[4.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT5_EM: _b[5.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT6_EM: _b[6.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT8_EM: _b[8.cat_chefe_agg] / _b[7.cat_chefe_agg]), post
estimates store nlcom_ens_medio_cat_2011

*2012
quietly logistic ens_medio b(7).cat_chefe_agg b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2012) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_cat_2012

quietly margins cat_chefe_agg, post
estimates store margins_ens_medio_cat_2012

quietly nlcom (CAT1_EM: _b[1.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT2_EM: _b[2.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT3_EM: _b[3.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT4_EM: _b[4.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT5_EM: _b[5.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT6_EM: _b[6.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT8_EM: _b[8.cat_chefe_agg] / _b[7.cat_chefe_agg]), post
estimates store nlcom_ens_medio_cat_2012

*2013
quietly logistic ens_medio b(7).cat_chefe_agg b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2013) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_cat_2013

quietly margins cat_chefe_agg, post
estimates store margins_ens_medio_cat_2013

quietly nlcom (CAT1_EM: _b[1.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT2_EM: _b[2.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT3_EM: _b[3.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT4_EM: _b[4.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT5_EM: _b[5.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT6_EM: _b[6.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT8_EM: _b[8.cat_chefe_agg] / _b[7.cat_chefe_agg]), post
estimates store nlcom_ens_medio_cat_2013

*2014
quietly logistic ens_medio b(7).cat_chefe_agg b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2014) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_cat_2014

quietly margins cat_chefe_agg, post
estimates store margins_ens_medio_cat_2014

quietly nlcom (CAT1_EM: _b[1.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT2_EM: _b[2.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT3_EM: _b[3.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT4_EM: _b[4.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT5_EM: _b[5.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT6_EM: _b[6.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT8_EM: _b[8.cat_chefe_agg] / _b[7.cat_chefe_agg]), post
estimates store nlcom_ens_medio_cat_2014

*2015
quietly logistic ens_medio b(7).cat_chefe_agg b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2015) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_cat_2015

quietly margins cat_chefe_agg, post
estimates store margins_ens_medio_cat_2015

quietly nlcom (CAT1_EM: _b[1.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT2_EM: _b[2.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT3_EM: _b[3.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT4_EM: _b[4.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT5_EM: _b[5.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT6_EM: _b[6.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT8_EM: _b[8.cat_chefe_agg] / _b[7.cat_chefe_agg]), post
estimates store nlcom_ens_medio_cat_2015

*2016
quietly logistic ens_medio b(7).cat_chefe_agg b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2016) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_cat_2016

quietly margins cat_chefe_agg, post
estimates store margins_ens_medio_cat_2016

quietly nlcom (CAT1_EM: _b[1.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT2_EM: _b[2.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT3_EM: _b[3.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT4_EM: _b[4.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT5_EM: _b[5.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT6_EM: _b[6.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT8_EM: _b[8.cat_chefe_agg] / _b[7.cat_chefe_agg]), post
estimates store nlcom_ens_medio_cat_2016

*2017
quietly logistic ens_medio b(7).cat_chefe_agg b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2017) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_cat_2017

quietly margins cat_chefe_agg, post
estimates store margins_ens_medio_cat_2017

quietly nlcom (CAT1_EM: _b[1.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT2_EM: _b[2.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT3_EM: _b[3.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT4_EM: _b[4.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT5_EM: _b[5.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT6_EM: _b[6.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT8_EM: _b[8.cat_chefe_agg] / _b[7.cat_chefe_agg]), post
estimates store nlcom_ens_medio_cat_2017

*2018
quietly logistic ens_medio b(7).cat_chefe_agg b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2018) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_cat_2018

quietly margins cat_chefe_agg, post
estimates store margins_ens_medio_cat_2018

quietly nlcom (CAT1_EM: _b[1.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT2_EM: _b[2.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT3_EM: _b[3.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT4_EM: _b[4.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT5_EM: _b[5.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT6_EM: _b[6.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT8_EM: _b[8.cat_chefe_agg] / _b[7.cat_chefe_agg]), post
estimates store nlcom_ens_medio_cat_2018

*2019
quietly logistic ens_medio b(7).cat_chefe_agg b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2019) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_cat_2019

quietly margins cat_chefe_agg, post
estimates store margins_ens_medio_cat_2019

quietly nlcom (CAT1_EM: _b[1.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT2_EM: _b[2.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT3_EM: _b[3.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT4_EM: _b[4.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT5_EM: _b[5.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT6_EM: _b[6.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT8_EM: _b[8.cat_chefe_agg] / _b[7.cat_chefe_agg]), post
estimates store nlcom_ens_medio_cat_2019

*2020
quietly logistic ens_medio b(7).cat_chefe_agg b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2020) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_cat_2020

quietly margins cat_chefe_agg, post
estimates store margins_ens_medio_cat_2020

quietly nlcom (CAT1_EM: _b[1.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT2_EM: _b[2.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT3_EM: _b[3.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT4_EM: _b[4.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT5_EM: _b[5.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT6_EM: _b[6.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT8_EM: _b[8.cat_chefe_agg] / _b[7.cat_chefe_agg]), post
estimates store nlcom_ens_medio_cat_2020

*2021
quietly logistic ens_medio b(7).cat_chefe_agg b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2021) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_cat_2021

quietly margins cat_chefe_agg, post
estimates store margins_ens_medio_cat_2021

quietly nlcom (CAT1_EM: _b[1.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT2_EM: _b[2.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT3_EM: _b[3.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT4_EM: _b[4.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT5_EM: _b[5.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT6_EM: _b[6.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT8_EM: _b[8.cat_chefe_agg] / _b[7.cat_chefe_agg]), post
estimates store nlcom_ens_medio_cat_2021

*2022
quietly logistic ens_medio b(7).cat_chefe_agg b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2022) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_cat_2022

quietly margins cat_chefe_agg, post
estimates store margins_ens_medio_cat_2022

quietly nlcom (CAT1_EM: _b[1.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT2_EM: _b[2.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT3_EM: _b[3.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT4_EM: _b[4.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT5_EM: _b[5.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT6_EM: _b[6.cat_chefe_agg] / _b[7.cat_chefe_agg]) (CAT8_EM: _b[8.cat_chefe_agg] / _b[7.cat_chefe_agg]), post
estimates store nlcom_ens_medio_cat_2022


****************Tabelas

*OR
etable, estimates(mod_ens_medio_cat_1992 mod_ens_medio_cat_1993 mod_ens_medio_cat_1995 mod_ens_medio_cat_1996 mod_ens_medio_cat_1997 mod_ens_medio_cat_1998 mod_ens_medio_cat_1999 mod_ens_medio_cat_2001 mod_ens_medio_cat_2002 mod_ens_medio_cat_2003 mod_ens_medio_cat_2004 mod_ens_medio_cat_2005 mod_ens_medio_cat_2006 mod_ens_medio_cat_2007 mod_ens_medio_cat_2008 mod_ens_medio_cat_2009 mod_ens_medio_cat_2011 mod_ens_medio_cat_2012 mod_ens_medio_cat_2013 mod_ens_medio_cat_2014 mod_ens_medio_cat_2015 mod_ens_medio_cat_2016 mod_ens_medio_cat_2017 mod_ens_medio_cat_2018 mod_ens_medio_cat_2019 mod_ens_medio_cat_2020 mod_ens_medio_cat_2021 mod_ens_medio_cat_2022) cstat(_r_b, nformat(%4.3f)) stars( , attach(_r_b) ) mstat(N, nformat(%9.0g))

collect layout (coleq#colname#result[_r_b] result[N] result[aic] result[bic] result[chi2] result[df_m] result[p] result[ll] result[r2_p]) (etable_depvar#stars) (), name(ETable)
collect export "C:\Users\andre\Desktop\Artigos\Artigo_Ens_Sup_2\Resultados\Modelos.xlsx", as(xlsx) sheet(Ens_Med_cat_OR, replace) cell(A1) modify


*RR
etable, estimates(nlcom_ens_medio_cat_1992 nlcom_ens_medio_cat_1993 nlcom_ens_medio_cat_1995 nlcom_ens_medio_cat_1996 nlcom_ens_medio_cat_1997 nlcom_ens_medio_cat_1998 nlcom_ens_medio_cat_1999 nlcom_ens_medio_cat_2001 nlcom_ens_medio_cat_2002 nlcom_ens_medio_cat_2003 nlcom_ens_medio_cat_2004 nlcom_ens_medio_cat_2005 nlcom_ens_medio_cat_2006 nlcom_ens_medio_cat_2007 nlcom_ens_medio_cat_2008 nlcom_ens_medio_cat_2009 nlcom_ens_medio_cat_2011 nlcom_ens_medio_cat_2012 nlcom_ens_medio_cat_2013 nlcom_ens_medio_cat_2014 nlcom_ens_medio_cat_2015 nlcom_ens_medio_cat_2016 nlcom_ens_medio_cat_2017 nlcom_ens_medio_cat_2018 nlcom_ens_medio_cat_2019 nlcom_ens_medio_cat_2020 nlcom_ens_medio_cat_2021 nlcom_ens_medio_cat_2022) cstat(_r_b, nformat(%4.3f)) cstat(_r_lb, nformat(%4.3f)) cstat(_r_ub, nformat(%4.3f)) export("C:\Users\andre\Desktop\Artigos\Artigo_Ens_Sup_2\Resultados\Modelos.xlsx", as(xlsx) sheet(Ens_Medio_cat_RR, replace) cell(A1) modify)
























